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I.  INTRODUCTION 


Projectile  impact  shock  initiation  of  high  explosive  has  long 
been  a  subject  of  considerable  interest  in  the  energetic 
materials  community.  Considerable  experimental  data  has  been 
generated  over  the  years.  Numerical  modeling  of  projectile 
impact  shock  initiation  for  comparison  with  experiments  has  been 
reported  in  at  least  one  case.^  However,  a  detailed  analysis  of 
the  flow  fields  revealed  by  the  computations  was  not  presented. 

Accordingly,  we  consider  in  this  work  the  problem  of 
oomposition-B  (bare  as  well  as  confined)  impacted  by  small 
hi gh-speedr steel  projectiles.  Using  the  reactive  hydrodynamic 
codes  S I II J  and  2  D  P.  we  have  addressed  the  posed  problem 
numerically.  This  approach  is  favored  by  virtue  of  its 
simulation  capability,  detailed  results,  moderate  cost,  and  lack 
of  risk  to  the  investigator.  In  the  following  section,  we 
describe  the  numerical  simulation  of  projectile  impact,  on  bare 
charges  for  six  different  projectile  diameters.  Results  obtained 
with  2DE  for  projectile  impact  arc  compared  with  planar  impact 
results  obtained  using  SIM.  The  results  were  used  to  generate 
linear  fits  which  suffice  to  describe  the  impact  shock  initiation 
of  bare  charges  consistent  with  the  critical  energy  criterion. 


1.  D.  C.  Slade  and  J .  Dewey,  "High-Order  Initiation  of  Two  Military  Explosives 
by  Progeotile  Impact,"  Ballistic  Research  Laboratory  Report  No.  1021 , 

July  1957  (AD  145868 ). 

2.  S.  M.  Brown  and  E.  G.  Whitbread,  " The  Initiation  of  Detonation  by  Shock 
Waves  of  Known  Duration  and  Intensity ,"  Les  Ondes  De  Detonation,  C.N.R.S. 

No.  109  (Paris,  1962),  pp  69-80. 

3.  L.  A.  Roslund,  J.  W.  Watt,  and  N.  L.  Coleburn,  " Initiation  of  Warhead 
Explosives  by  the  Impact  of  Controlled  Fragments  I.  Normal  Impact,  " 

Naval  Ordnance  Laboratory  Technical  Report  NOLTR-7 3-124,  August  1974. 

4.  K.  L.  Bahl,  H.  C.  Vantine  and  R .  C.  Wllingort,  " The  Shock  Initiation  of 
Bare  and  Covered  Explosives  by  Projectile  Impact Seventh  Symposium 
(International)  on  Detonation,  June  1981 ,  pp  325-335. 

5.  C .  L.  Mader  and  M.  S.  Shaw ,  " Users  Manual  for  SIN,  A  One-Dimensional 
Hydrodynamic  Code  for  Problems  which  Include  Chemical  Reactions,  Elastic 
Plastic  Flow,  Spalling  Phase  Transitions,  Melting,  Forest  Fire,  Detonation 
Buildup,  and  SESAME  Tabular  Equation  of  State,  "  Los  Alamos  Scientific 
Laboratory  Report  LA-7264-M,  September  1978. 

6.  J .  D.  Kershner  and  C.  L.  Mader,  tf2DE,  A  Two-Dimensional  Continuous  Eulerian 
Hydrodynamic  Code  for  Computing  Multicomponent  Reactive  Hydrodynamic 
Problems,"  Los  Alamos  Scientific  Laboratory  Report  LA-4846,  March  1972. 
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We  then  turned  our  attention  to  a  numerical  assessment  of  shock 
initiation  in  covered  charges.  The  results  show  the  buildup  of 
shock  to  detonation  in  detail  by  means  of  contour  p^ots.  We  also 
checked  our  results  with  the  Jacobs-Roslund  formula"' 


II.  PROJECTILE  IMPACT  SHOCK  INITIATION  OF  BARE  EXPLOSIVE  CHARGES 
Geometry  and  Computational  Considerations 

In  our  computations,  wc  have  considered  steel  cylinders  of' 
aspect  ratio  (£,/$)  equal  to  one  and  corn  position-  B  charges  of 
aspect  ratio  one- half.  Sufficient  target  material  is  provided 

when  the  explosive  diameter  is  three  tines  the  projectile 
diameter.  This  aspect  ratio  is  considered  adequate  to  take 

account  of  rarefaction  effects  from  free  surfaces  and  edges. 
Moreover,  Brown  and  Whitbread^  show  that  different  aspect  ratios, 
except  for  the  case  of  a  disc  ({,/ <  1/4  )  ,  do  not  result  in 
different  critical  velocities  for  shock  initiation. 

With  the  considerations  mentioned  above,  we  have  set  up 
impact  problems  for  2DE  calculations  with  axi  symmetric  grids  as 
summarized  in  Table  1.  Here  <J>  is  the  projectile  diameter,  AR 
the  radial  cell  size,  I  the  number  of  cells  along  the  radial 
axis,  A  Z  the  axial  cell  size,  J  the  number  of  cells  along  the 
axis  of  symmetry,  AT  the  time  step  of  each  computing  cycle,  and  N 
the  total  number  of  cycles  to  be  completed. 

TABLE  1 


( tom ) 

Input 

AR 

(mm) 

Data  for 

I 

2DE  Computational 

AZ 

( mm ) 

Grids 

J 

-  Hare  Charges 

AT 

(  ys) 

N 

5 

0.25 

45 

0 . 25 

6  8 

0 . 005 

4  00 

3 

0.40 

4  5 

0.40 

63 

0.010 

250 

1  n 

0 .50 

4  5 

0.50 

6  8 

0.010 

4  00 

1  2 

Q  .  4  0 

65 

0.40 

98 

0.010 

530 

1  5 

0 .50 

6  5 

0 .50 

9  8 

0.010 

530 

1  8 

0.60 

C  5 

0.60 

9  8 

0.015 

4  00 
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Results 


! 

A  number  of  graphical  representations  of  our  numerical 
results  are  available.  The  sequence  of  events  in  projectile 
impact  shock  initiation  |s  most  clearly  illustrated  in  the 
series  of  mass  fraction  contour  plots  of 'Figure  la.  The  mass 
fraction  varies  from  one  to  zero  as  chemical  reaction  in  the 
explosive  runs  to  completion.  The  plots  show  results  for  impact 
of  a  5  mm  diameter  projectile  at  1.4  mrn/ys.  The  corresponding 
isobar  plots  are  s h o w n  in  Figure  1b.  Detonation,  which'  may  be 
recognized  by  the  close  spacing  of  the  contour  lines,  is  observed 
to  begin  at  the  shock  front  and  spread  outward .  Similar  plots 
for  the  5  mm  projectile  at  a  subcritical  impact  velocity  of  1.2 
rrn/ys  are  shown  in  Figures  2a  and  2b.  In  this  case,  no 
detonation  develops  and  concentration  of  the  mass  fraction 
contour  lines  is  not  observed.  Figure  3  shows  mass  fraction 
contour  plots  for  a  projectile  15  mm  in  diameter  with  an  impact 
velocity  of  0.7  mm/jus,  in  which  case  initiation  resulted.  In 
contrast  to  the  small  diameter,  high  velocity  case,  initiation 
here  first  appears  after  the  shock  has  propagated  some  distance 
from  the  impact  point.  The  ensuing  detonation  propagates 
outwards.  This  latter  observation  is  typical  of  all  the 
computations  for  Kj)  >  5  mm. 


We  also  made  computations  comparing  projectiles  of  different 
aspect  ratio.  Mass  fraction  contour  plots  for  the  1.2  inm/ys 
impact  of  a  10  mm  diameter  projectile  with  aspect  ratios  of  1.0 
and  1.5  are  shown  in  Figure  4.  The  results  show  that  the  flow 
fields  created  in  the  two  cases  are  essentially  identical  during 
the  period  required  for  buildup  to  detonation.  Computations  were 
made  with  impact  velocity  varied  in  steps  of  0,2  rnm/ys  until  the 
minimum  velocity  for  a  shock-to-detonation  transition  was 
determined  for  each  projectile  diameter.  Some  one-  and 

two-dimensional  results  are  listed  in  Table  2,  where  $  is  the 
projectile  diameter  in  the  case  of  two-dimensional  computations 
and  flyer  plate  thickness  in  the  case  of  one-dimensional 
computations,  V  *  is  the  critical  velocity,  P  ^  the  peak  pressure 
of  the  initial  shock  pulse  in  the  explosive,  x#  the  distance  of- 
shock  run  to  detonation,  T#  the  time  of  run  to  detonatiion,  and 
is  a  parameter  equivalent  to  the  critical  energy  criterion.  We 
also  include  the  corresponding  results  p.,  x  *  ,  and  t*.  obtained 
from  the  one-dimensional  computation  using  SIN  which  produces  the 
same  critical  velocities.  Using  the  detonation  pressure,  PpT  = 
23.4  GPa,  for  composition-B,  we  have  calculated  the  average  rate 
of  shock  buildup  Y>)t,  y^  and  the  gradient  of  shock  buildup  Z ^  ,  z* 
defined  as  follows: 


x  * 
y* 


=  (p 
=  (p 
=  (p 
=  (p 


C  J 
C  J 
C  J 


-  Pt)/T*, 

-  Pp/t*, 

-  pp/x*’ 

-  p  r )  /  x , 


It  should  be  noted  that1conditions  cited 
nay  be  very  close  to  critical  or  as  much 
cited  value.  Thus,  our  definition  of 
entirely  uniform. 


as  critical  in  the  table 
as  0.2  mm/ys  above  the 
critical  values  is  not 
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Figure  la.  Mass-Fraction  Contours  Showing  the  Spread  of  Reaction  in 
Bare  Charge  ((J>  =  5  mm,  V  *  1,4  mm/ys  )  . 
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Figure  lb.  Isobars  Showing  Shock  and  Reaction  in  Bare  Charge 
(4>  =  5  mm,  V  =  1.4  mm/ys) . 
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Figure  2a. 


Mass-Fraction  Contours  for  "NO  GO"  in  Bare  Charge 
(<j>  =  5  ram,  V  =  1.2  mm/ys) . 
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Figure  2b.  Isobars  for  "NO  GO"  in  Bare  Charge 
(<f)  =  5  mm,  V  =  1.2  mm/ys) . 
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Figure  3.  Mass-Fraction  Contours  Showing  the  Spread  of  Reaction  in 
Bare  Charge  (()>  =  15  nun,  V  =  0.7  mm/ys)  . 
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Figure  31.  Mass-Fraction  Contours  (continued) 
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Figure  4.  Comparison  of  Different  Aspect  Ratios  (cf> 
V  =  1.2  mm/ys) . 


Summary  of  Computed  Critical  Values  -  Bare  Charges. 
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In  Figure  5  we  have  plotted  'critical  velocity  obtained  in  our 
2DE  calculations  versus  the  parameter  l//<j>  for  comparison  with 
the  NOL  empirical  formula.  Our  computed  critical  velocities  fall 
slightly  below  the  experimental  values  but  still  exhibit  the 
characteristic  linear  dependence  on  i//<j>  .  The  straight  line  of 
Figure  5  can  be  described  by  the  equation: 

;  i 

\V*  =  ai  ;+  bl  //<()  CD 

3/2 

with  constants  a  1  =  -0.16  mm/ys  and  b  =  3-33  mm  /ys.  The 
critical  velocity  at  cj>  =  5mm  appears  to  deviate  somewhat  from 
the  linear  relation  and  has  not  been  included  in  the  fit.  Thi$ 
is  apparently  a  consequence  of  the  alternate  mode  of  buildup 
observed  for  the  small— diameter  high-velocity  case.  It  in  ay  be 
noted  that  Equation  (1)  has  been  used  for  a  long  time.  (Note  a  =o 
and  b^=3.26  mil  2  /y  s  for  composition-B  in  Reference  1.)  ' 

As  previously  noted,  our  cited  critical  conditions  vary  from 
the  actual  conditions  in  a  nonuniform  manner.  It  is  of  interest 
to  observe  how  closely  the  paths  to  initiation  conform  to  the 
Pop-plot  for  comp-B.  This  is  illustrated  in  Figure  6.  The  SIN 
results  are  all  seen  to  lie  reasonably  close  to  the  Pop-plot 
indicating  that  the  single-curve  buildup  hypothesis  (implied  in 
Forest  Fire)  is  applicable  to  the  one-dimensional  situation. 
However,  for  the  multidimensional  projectile  impact  case  buildup 
deviates  considerably  from  the  Pop-plot.  This  appears  to  be 
chiefly  due  to  the  near  critical  nature  of  the  cases  considered 
while  the  Pop-plot  represents  supercritical  buildup.  For 
example,  the  computation  at  ^  =  10  mm  almost  doesn't  initiate. 
The  pressure  first  builds,  then  drops  off  and  finally  rises  again 
to  the  CJ  value.  This  is  clearly  our  most  nearly  critical  case. 
The  computation  atcj>  =  8  mra  i  s  more  supercritical  and  conforms 
more  exactly  to  the  Pop-plot.  Only  one  point  could  be  obtained 
at<f>  =  10  mm.  This  point  lies  below  the  Pop-plot  while  the 
results  at  larger  diameters  were  generally  above  the  Pop-plot. 
This  gives  a  third  indication  that  the  buildup  to  detonation 
proceeds  somewhat  differently  in  the  small-diameter, 
high-velocity  case. 

In  the  literature,  the  critical  energy  "Ctiterion  ^  has  been 
very  popular  for  the  study  of  one-dimensional  shock  initiation 
(as  in  the  wedge  test).  If  this  criterion  is  applicable  to  our 
two-dimensional  study,  then  we  must  have  V*  <p  =  constant.  Yet  our 
numerical  results  of  Table  2  do  not  maintain  a  strict  constancy 
°f  V*  <j>  .  However,  some  of  this  may  be  explained  by  the 
nonuniform  determination  of  the  critical  condition  . 

In  the  original  derivation  of  the  critical  energy  criterion, 
the  shape  and  width  of  the  initial  shock  pulse  play  an  important 
role  in  characterizing  the  shock  initiation.  In  our 

7.-  F.  E.  Walker  and  R.  J.  Wasley,  "  Critical  Energy  for  Shock  Initiation  of 
Heterogeneous  Explosives,  "  Explosivstoffe  Mr.  1,  1969,  pp  9-13, 
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CRITICAL  VELOCITY  V.  | 


Figure  5.  Critical  Velocity  vs  l//§  for  Bare  Explosive  Charges. 
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Figure  6.  Comparison  of  Shock  Buildup  Trajectories  to  Pop-Plot- 
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two-dimensional  Impact  problem,  the  pulse  shape  and  duration  jr? 
jffooted  by  Tftgefftctioti  waves  issuing  from  free  surfaces  and 
edges.  It  appears,  however,  that  the  reactive  shock  in  the 
explosive  can  accelerate  to  detonation  by  self-support  provided 
the  necessary  conditions  are  satisfied.  A  number  of  other  linear 
relations  among  the  problem  parameters  were  observed  t.  ;i 
illustrated  in  Figures  7-11.  The  following  linear  fits  can  S"rvc 
to  complement  Equation  ( 1 )  and  provide  n  full  description  of  the 
initiation  process  in  terms  of  the  stimulus. 


j*  ^ 

Pj-  -  +  b2V*  =  9.2  +  (2) 

T*  =  a3  +  b3^  (3) 

X*  =  a4  +  b4T*  (4) 

>  > 

Y*  =  a5  _+  b5V*  =  a5  +  b5//<j)  (5) 

z*  *  a6  +  b6V*  =  a6  +  b6//<j>  (6) 


The  values  of  the  constants  are  given  in  Table  3.  It  should  be 
noted  that  Equation  (4)  represents  a  relationship  that  is  a 
property  of  the  explosive  only.  The  foregoing  linear  relations 
also  appear  to  hold  good  for  our  SIN  results.  In  fact,  the  SIN 
results  fit  straight  lines  even  better  in  Figures  7-11. 


III.  SHOCK  INITIATION  CF  COVERED  EXPLOSIVE  CHARGES 

We  have  also  addressed  the  related  problem  posed  by 
introducing  a  steel  plate  of  thickness  h  between  the  steel 
projectile  and  the  composition-B.  Table  4  summarizes  our  new 
grids.  We  have  considered  projectile  diameters  of  5  ♦  6  ,  1 0 f 1 2  and 
15  mm  and  cover  plates  of  one-third  and  one-half  the  diameter  in 
each  case. 


Detonation  in  the  covered  charges  develops  in  a  manner 
similar  to  the  bare  charges  as  illustrated  in  the  mass-fraction 
contour  plots  of  Figure  12  for  a  15  tn in  projectile  with  an  impact 
velocity  of  1.0  mm /pis. 

The  critical  pa r ame t e r - v a lu e s  are  summarized  in  Table  5.  For 
h/  4>  =  1/3,  the  critical  energy  criterion  ( y2  ^  _  Const)is  followed 
almost  exactly  except  in  the  5  mrn  diameter  *case.  The  constancy 

of  V*(f>,  is  also  good  for  h/d>  =  1/2. 
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Figure  7.  Linear  Fits  of  Initial  Shock  Pressure. 


Figure  8. 


Linear  Fit  of  Time  to  Detonation  versus  Projectile 
Diameter, 


Figure  9.  Linear  Fit  of  Distance  to  Detonation  versus  Time  to 
Detonation . 
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TABLE  4 


Input  Data  For  2  D  E  Computational  Grids 


Covered  pharges 


h/  <(> 

♦ 

(mm) 

AR 

(mm) 

I 

AZ 

(mm) 

J 

AT 

(PS) 

N 

5 

0.125 

45 

0.125 

125 

0 .0025 

1  ,  000 

8 

0.200 

4  5 

0 .200 

■  125 

0.0050 

90  0 

1  /3 

1  0 

0.250 

45 

0.250 

125 

0 .0050 

0  00 

1  2 

0.240 

55 

0.240 

1  40 

0 .0050 

900 

1  5 

0  .  300 

55 

0 . 300 

1  4  0 

0 . 0020 

700 

5 

0.125 

45 

0.125 

1  4  0 

0 .0025 

1  ,  000 

8 

0 . 200 

4  5 

0 .200 

125 

0 .0050 

1  ,  ooo 

1  /2 

1  0 

0.250 

45 

0.250 

125 

0 .0050 

8  0  0 

12 

0.240 

55 

0.240 

1  40 

0.0050 

1,000 

15 

0.300 

55 

0 . 300 

1  4  0 

C  .  0  0  P  0 

700 

h/* 

Summary  of 

♦  V 
(mm)  (mm/ps) 

TABLE  5 

Computed  Critical  Values  - 

V*<f>  Pj  x* 

(mm'Vps2)  (GPa)  (mm) 

Covered 

T* 

(ps) 

Char  g  e  3 

Y*  Z* 

(GPa/ps) (GPa/ mm) 

5 

1  .  4 

9.80 

10.0 

1.13 

0 .24 

76.6 

16.3 

8 

1  .  2 

12.10 

10.0 

1  .  8 

0.35 

5  2.6 

10.2 

1/3 

1  0 

1  .  1 

12.10 

10.0 

2.3 

0.50 

36.8 

8 . 0 

1  2 

1  .  0 

12.00 

8 . 0 

2.5 

0.00 

34.0 

£  .  1 

15 

0.9 

12.15 

7.0 

3.6 

0.90 

2  3.8 

5.9 

5 

1.7 

14.45 

14.0 

0.63 

0. 1 

14  4.0 

22.9 

8 

1  .  4 

15.60 

14.0 

0. 8 

0.3 

48.0 

18.0 

1/2  10 

1  .  3 

15.90 

12.0 

1  .  0 

0.3 

54.7 

16.4 

1  2 

1  .  2 

16.15 

11.0 

2.3 

0. 4 

4  3.5 

7.6 

15 

1  . 0 

15.00 

10.0 

3.6 

1  .  0 

18.4 

5.  1 
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Mass-Fraction  Contours  Showing 
Covered  Charge  (4)  =  15  mm,  h/<J> 
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The  Jacobs-Roslund  ( N  0  L )  formula'  for  normal  impact  against 
eovered  eharges  may  be  written 

V*  =  V*  (1  +  h/ <J>) _1  =  Vo  +  A/4  (7) 

>* 

Plotting  V*  versus  l//<(>  ,  as  in  Figure  13,  allows  a  comparison  of 
the  computed  results  ath/ij)  =  0,l/3and  1/2  with  the  results  of  the 
Jaeobs-Roslund  formulas  with  V  =0  and  A  =  3.13  mm  .  Only 

at  the  smaller  projeetile  diameSer  do  the  results  for  the  covered 
charges  differ  significantly  from  the  linear  relationships 
defined  by  the  bare  eharge  results  .  In  addition,  the  covered 
charge  results  tend  to  diverge  from  one  another  at  smaller 
diameters  . 


Plots  corresponding  to  those  of  Figures  7  —  11  for  the  bare 
eharge  problem  have  also  been  constructed  in  the  eovered  eharge 
case.  These  are  shown  in  Figures  14-17.  Figure  16,  which 
corresponds  to  Figure  Q,  shows  the  distance  to  detonation  plotted 
versus  the  time  to  detonation.  As  previously  noted,  this  is 
strictly  an  explosive  property  oeeording  to  the  single-curve 
buildup  hypothesis  and  should  be  independent  ofh/<f>. 


IV.  SUMMARY 

Using  the  SIN  and  2DE  hydrodynamic  computer  eodes,  vie  have 
conducted  a  numerical  study  of  the  shock  initiation  response  of 
both  bare  and  covered  eonposi t i on-B  eharges  to  projeetile  impact. 

In  the  bare  charge  ease  we  observed  two  modes  of  shoek 
initiation.  In  the  ease  of  larger  projeetile  diameters,  the 
detonation  is  observed  to  begin  at  some  distance  from  the  impact 
point  and  propagate  outward.  When  the  projeetile  diameter  is 
sufficiently  small  (and  the  eritieal  velocity  correspondingly 
large)  detonation  oecurs  almost  immediately.  We  also  verified 
t  hi  e  independence  of  the  shoek  initiation  response  with  respeet  to 
projeetile  aspeet  ratio.  We  then  determined  eritieal  conditions 
for  shock  initiation  for  six  different  projectile  diameters  with 
unit  aspect  ratio.  Except  for  the  case  of  the  smallest 

projectile  ($=  5  mm),  the  eritieal  velocity  was  observed  to  have 
a  linear  dependence  on  l//cf)  with  values  falling  a  little  below 
the  NOL  empirical  line.  The  progress  of  buildup  to  detonation 
was  observed  and  found  to  elosely  folio w  the  Pop-plot  when  the 

impact  was  sufficiently  supereritieal.  Finally,  we  generated 
additional  linear  relations  to  completely  describe  buildup  to 
detonation  in  the  eritieal  ease. 

In  the  ease  of  eovered  charges,  we  found  that  the  development 
of  detonation  is  similar  to  that  observed  for  bare  eharges, 
Critieal  velocities  and  other  parameters  were  determined  as  for 

the  bare  charges.  The  critical  velocities,  when  adjusted  for 
cover  plate  thickness  according  to  the  Jaeobs-Roslund  formula, 
were  found  to  agree  well  with  the  bare  eharge  results  at  larger 
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Figure  13.  Comparison  of  Bare  and  Covered  Charge  Critical  Velocity  with  Jacobs-Roslund 
Formula, 
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Figure  14.  Linear  Fits  of  Initial  Shock  Pressure. 


Figure  15.  Linear  Fit  of  Time  to  Detonation  versus  Projectile  Diameter. 


Figure  16.  Linear  Fit  of  Distance  to  Detonation  versus  Time  to  Detonation. 
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Figure  17.  Linear  Fits  of  Shock  Buildup  Gradient. 


projectile  diameters  and  to  diverge  from  the  bare  charge  results  and  from  each 
other  with  decreasing  projectile  diameter. 
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